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ABSTRACT 

■ We present three dimensional smoothed particle hydrodynamics (SPH) calculations 

of irradiation-driven warping of accretion discs. Initially unwarped planar discs are 
unstable to the radiation reaction when the disc is illuminated by a central radiation 
source. The disc warps and tilts and precesses slowly in a retrograde direction; its 
shape continuously flexes in response to the changing orientation of the Roche poten- 
tial. We simulate ten systems: eight X-ray binaries, one cataclysmic variable (CV), 
Jjj ' and a 'generic' low mass X-ray binary (LMXB). We adopt system parameters from 

observations and tune a single parameter: our model X-ray luminosity (£*) to repro- 
duce the observed or inferred super-orbital periods. Without exception, across a wide 
range of parameter space, we find an astonishingly good match between the observed 
Lx and the model . We conclude irradiation-driven warping is the mechanism un- 
derlying the long periods in X-ray binaries. Our HerX-1 simulation simultaneously 
| reproduces the observed Lx , the "main-" and "short-high" X-ray states and the 

04 ■ orbital inclination. Our simulations of SS433 give a maximum warp angle of 18.6°, 

' a good match to the cone traced by the jets, but this angle is reached only in the 

outer disc. In all cases, the overall disc tilt is less than 13° and the maximum disc 
warp is less than and or equal to 21°. In particular, the disc warp in 4U 1626-67 can- 
not explain the observed torque reversals. Taking our results at face value, ignoring 
<^> the finite opening angle of the disc, we deduce orbital inclinations of approximately 

77° for 4U 1916-053 and approximately 69° for 4U 1626-67. We also simulate CygX-2, 
SMCX-1, CygX-1, and LMCX-3. For high mass X-ray binary (HMXB) parameters, 
the discs' maximum angular elevation is invariably at the outer edge. For LMXBs 
' with extreme mass ratios a strong inner disc warp develops, completely shadowing 

\ parts of the outer disc. This inner warped disc executes retrograde precession while 

the outer disc executes prograde apsidal precession. The remaining LMXBs develop 
a less extreme warp in the inner disc, with the entire disc tilting and precessing in 
a retrograde direction. For our CV, KRAur, we matched the inferred disc precession 
period by adopting Lx — 10 37 erg sec -1 , which would require steady nuclear burning 
on the white dwarf surface. 

Key words: accretion: accretion discs - X-rays: binaries - binaries: close - X-rays: 
binaries - methods: numerical - Stars: individual: Her X-l 



1 INTRODUCTION 

In iFoulkes et al. I (|2006l . hereafter Paper I), we presented 
three dimensional smoothed particle hydrodynamics (SPH) 
simulations of warped accretion discs in X-ray binary sys- 
tems. We showed that geometrically thin, optically thick, 
initially planar, accretion discs when irradiated by a cen- 
tral radiation source are unstable to warping. The central 
illumination of the non-axisymmetric disc surface causes a 
radiation pressure torque that acts on the disc to induce a 



twist or warp. In general the warps precess in a retrograde 
direction. Paper I discussed X-ray binary systems with a 
range of mass ratios and luminosities. Here we extend the 
work of Paper I, and simulate the well-known binary systems 
Hercules X-l (Her X-l), SS433, several other X-ray binary 
systems, and one cataclysmic variable (CV). The X-ray bi- 
naries modelled all have a super-orbital period, i.e. a period 
in the light curve that is much longer than the orbital period. 
We found for disc warping to occur the radiation intensity 
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Table 1. Major parameters for binary system HZHER/HerX-1 



Parameter 


value 




reference 


orbital period 


-^orb " 


= 1.7 days 


1 


binary separation 


a = 6 


■ 10 n cm 


2 


orbit eccentricity 


e = 3 


■ 10" 4 


2 


inclination of system 


i = 8( 


)° - 90° 


2,3,4 


mass of neutron star 


Mi = 


1.3M Q 


2,3 


mass of donor 


M 2 = 


2.2M 


2.3 


mass transfer rate 


-M2 


= 10 17 ' 9 g sec" 1 


5 


X-ray luminosity 


L = 10 373 erg sec 1 


6 


super-orbital period 


^"prcc 


= 34.88 ±0.12 days 


9 


disc radius 


^°disc " 


= 1 - 2 ■ 10 11 cm 


10,11,12 



References: 1 Tananbaum et al. 1 Jl972h. 2 iDeeter et al. 




1981 


)■ 

), 


3 Middleditch & Nelson 1 Jl976h, 4 iGerend & B 


oynton 


1976 


5 Schandl 1 Jl996h. 6 iMcCrav et al. 1 dl982h, 7 


Triicmper et al. 


1978 


), 8 |Pravdo et al. 1 19771). 9 Giacconi et al. 


Jl973h, 10 IOke 


(1976 


), n ICrosa & Bovnton ( 1980), 12 Middleditch (1983) 



has t o be greater than a critical value fc.f. lWiiers Sz Pringle I 

1999); for low radiation levels warping does not oc cur. 

We use a non- linear SPH code based on that of lMurrav"! 
l| 19961 . Il998l ) but completely rewritten and extended as de- 
scribed in Paper I. Section [2] describes the two best-known 
systems; in sections [3] U and [S] we present our formulation 
of the problem using SPH; section [6] presents and discusses 
our results; and section [7] is devoted to conclusions. 



2 SYSTEMS MODELLED 

In total we present models of ten systems in this paper. 
Here in section [5] we discuss in detail HerX-1 and SS433, 
two archetypal warped disc X-ray binaries. HerX-1 is bright 
and well-constrained observationally, so it is important for 
assessment of our simulations. SS433 has been rather less 
definitively characterised, and has kinematic elements be- 
yond the scope of our model: namely the relativistic jets 
and the accretion disc wind. Nontheless the model is an in- 
teresting counterpoint to that of HerX-1. 

2.1 HerX-1 

HerX-1 was discovered by iTananbaum et al. I (| 19721 ). It is 
an intermediate mass X-ray binary consisting of an X-ray 
pulsar and a 13th m agnitude blue variable star HZ Her 
l|Davidsen et al. I Il972h . Sinusoidal variations and regular 
eclipses of the X-ray source indicates the system is viewed 
at a high inclination with an orbital period of 1.70 days. 
Assumin g HZ Her fills its Roche l obe and using the eclipse 
duration lTananbaum et al. I (|l972l ) estimated the mass ratio 
of the system to be Mdonor/M comp act ~ 1.9. In Table Q] we 
sum marise the paramet e rs we have adopted for HerX-1. 

ITananbaum et al. I l| 19721 ) and iGiacconi et al. I (Il973l ) 
reported a super-orbital X-ray cycle due to periodic ob- 
scurations of the X-ray source by the disc material. The 
34.88 ± 0.12 day cycle is comprised of on and off states 
lasting for 11-12 and 24 days respectively. I Jones fc Forman I 
(1 19761 ) observed during the 24 day off state and discovered 



the "short-high" emission which is approximately 30% of the 
pea k intensity of the "mai n-high" . 

IGiacconi et al. I (|l973l ) also reported regular dips in the 
X-ray light curve, the duration of which was about 1 — 3 
hours. The frequency of the dips was related to the sum of 
frequencies of the system orbital and the super-orbital fre- 
quency. They reported two types of dips. "Preeclipse" dips 
which appear just before an eclipse, at orbital phase (4>orb) 
approximately 0.8, which propagate back in super-orbital 
phase as the 35 day phase increases. These dips recur on a 
period of approximately 1.65 days. The second set of dips 
are known as the "anomalous dips". These dips occur at 
4>orb ~ 0.5 and are observed during the "short-high". The 
X-ray dips can be explained by a number of different mod- 
els, but they are certainly produced by occultation of the 
central X-ray source by the accreting disc material. 

It is generally accepted that the 35 da y super-orbita l 
perio d is a conseq uence of a twist ed/warped (|Kumar l ll986l . 
Il99d ) and tilted l|Roberts I ll974T ) accretion disc that pre- 
cesse s in a retrograde direction about the central neutron 
star (|Petterson II 19751 . 1 19771 . 1 19781 ). An observer sees the disc 
nearly edge on and as regions of the disc move above the 
orbital plane by more than th e disc inclina tion the central 
X-ray source will be occulted. iKatz I |l973) proposed that 
the tidal force of HZ Her acts on the tilted accretion disc 
and causes it to precess, and lLarwood et al. I (| 19961 ) find 
this mechanism produces the correct precession rate. Nei- 
ther of these studies, however, explain the origin of the tilt 
and twist of the disc. 

ISchandl fc Meyer I l| 19941 ) and ISchandl I (|l996l ) used a 
coronal disc wind model to generate a warped disc through 
asymmetric wind induced force on the disc surface. They 
modelled an optically thick, geometrically thin accretion disc 
that was heated by X-rays originating from the neutron star, 
which resulted in a hot corona above the disc. Near the disc 
the corona was hydrostatically layered and optically thick. 
However, in the outer regions the gas exceeds the escape ve- 
locity and leaves the gravitation potential of the binary sys- 
tem. This coronal wind generates a reaction force on the disc 
surface, which results in the outer disc regions warping and 
the inner parts being tilted but unwarped. For an orbital in- 
clination of 80° they were able to reproduce the "main-high" 
and "short-high". They also used a simple model to gener- 
ate the observed X-ray dips including the back propagation 
of the dips as the super-orbital precession proceeds. 

2.2 SS433 

SS433 is an unusual, well-studie d binary syst em a ssociated 
with t he su pernova re mnant W50 lClark I (|l984T ) ; see lMargon I 
l|l984h and iFabrika I l|2006h for reviews. Intensive observa- 
tions led to a model that includes an evolved binary with 
a large mass transfer rate and P OT b approximately 13 days. 
A Roche lobe-filling donor star loses mas s to a black hole 
l|Blundell. Bowler, fc Schmidtobreick]|2008l ) via an accretion 
disc. Some of the inflow is directed through the interaction 
of the disc and the compact object into oppositely-directed 
relativistic precessing jets. 

"Moving" emission lines have extreme radial velocities 
with a period of approximately 162 days. They are described 
by the "kinematical model" of the jets driven by a pre- 
cessing accretion disc (jFabian fc Rees ll 19791 : lMilgrorn~1ll979l : 
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Table 2. Major parameters adopted for binary system SS433 



Parameter 



value 



reference 



orbital period 
inclination of system 
system mass ratio 
mass of compact star 
mass of donor 
mass transfer rate 
X-ray luminosity 
super-orbital period 
disc nodding 
disc tilt 



•Porb = 13.1 days 
i = 78.83° 
q = 2 

Mi = 10 M e 
M 2 = 20 M© 
-M 2 = 10 



10 



,36 , 



M yr" 
1 



erg sec 
Pprec = 162.5 days 
P nod = 6.2877 days 
disc tilt « 20° 



1 
1 

3t 

3 

3 

6 

2, 

1 

1 

1 



4, 5 



t We note the recent paper Blunde ll. Bowler. fc Schm idtobrcick 
(|2008h favours Mi = 16M , M 2 = 22M Q , and consequently q = 
1.38. The simulation re ported here predated t his finding. 
References: 



' Brinkmann, Kawai. Matsuoka & Fink 1 


; Wiiers & Pringle 1 il999h. 


Kotani et al. 1 


J Marshall et al. 1 J2002J), IShklovskii 1 


19811). 



1984) 



1991) 



1996) 



lAbell fc Margonlll979T l. The "moving" lines are a lso seen in 
the Chandra X-ray spectra (jMarshall et al. 1120021 ). Another 
majo r periodicity foun d is the so-called "nodding" of the 
disc (|Katz et al. II1982T ). This is a 6.28 day period and is at- 
tributed to the donor star's passage through the nodal line 
of the titled disc which causes tidal deformation of the disc. 
In addition to the jets there is evidence of a very powerful 
outflow fro m the disc with a velocity of appr oximately 1500 
km s -1 (see Bcgclma n. King, fc Pringle 120061 . and references 
therein). Figure la of lBarnes et al. I (|20061 ) shows the major 
elements of SS433, which include the compact object, the 
Roche filling donor star, the accretion disc, the precessing 
relativistic jets and the outflow in a form of extend disc or 
disc wind. 

The precessional and orbital periods have provided a 
very s atisfactory model of the m otion of the jets dKatz et al. I 
1 19821 ; [Collins fc Garasi I ll994T l and the photometric varia- 
tions i Goranskii et al.Tl998l ). However, there are still many 
unanswered questions about the system. These include, the 
formati on process for the jets a nd the disc precession mech- 
anism (|Wiiers fc Pringle l(l999T l. In Table we give a sum- 
mary of the main parameters we adopted for SS 433. 



is the material derivative operator l|Batchelor I [l970) . p, v 
and P are the fluid density, velocity and pressure values 
respectively. The gravitational potential is "I>, the viscous 
force is f v and the force due to the radiation on the fluid is 
firr- The equation for the rate of change of thermal energy 
per unit mass can be expressed as 



4 NUMERICAL METHOD 



(4) 



We solve the basic equations in sect ion [3] using a SPH (jLucv I 
1 19771 ; iGingold fc Monaghan~lll977l ) computer code first de- 
veloped by Murray (1996, 1998). 

4.1 SPH formulation 

SPH uses smoothed quantities defined by the equation 



(/«>= / f(r')W(\r-r'\,h)dr', 



(5) 



here /(r) is an arbitrary function, W is a smoothing kernel, 
h is the smoothing length and the integration is over the 
entire space. A Monte-Carlo representation of ((5]l is used 
for each particle position 

JV 

</(*)) = V^/(rj)W(|ri- rj|, ft), (6) 
f-f Pi 

3=1 

where rrij is the mass of particle j, pj is the density of the 
particle j and the sum is over all particles N. We use a 
variable smoothing length such that each particle, i, has an 
associated smoothing length hi. The kernel is symmetrized 
with respect to particle pairs in order to conserve momen- 
tum. 



4-1.1 Equations of motion 

The equation of motion for each particle i is given by 

mi—rr = Fpa + Fa.i + F„; sc a + Fi rr A, (7) 
at 

which gives the force on particle i as a summation of the 
pressure, gravity, viscous and irradiation forces, Fp,i, Fg,„ 
F V i SC ,i and F irr ,i respectively. 



3 BASIC EQUATIONS 

We take the basic hydrodynamics equations governing the 
viscous gaseous material of the accretion disc. It is assumed 
that the disc self-gravity is negligible in comparison to the 
binary component's tidal field. The equations of continuity 
and motion may be written as 



Dt p 



V$ + f v + f ir 



where in the above 

D d „ 
— =. h v • V 

Dt dt 



(1) 
(2) 

(3) 



4-1.2 The pressure term 

For the pressure gradient we make use of the fact that 

1 ( P\ P 

-VP = V - +-Vp, 
P \P j P 2 



P J i ^ ' \Pi Pj 
Similarly the SPH energy equation is 



dui 1 



Pi . P , 



(8) 



and the pressure is given by 

-VP) =ym j (^ + ^) VW(|r,-rj|,ft). (9) 



,11 2 E m 4y + ^ Ivu-Vil-ln-rj./M. (10. 

3 3 
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4-1.3 The artificial viscosity term 

To stabilise the shock flo w regions we use d an artificial vis- 
cosity term as detailed by iMurravl l| 19961 ) and lTruss et al. I 
(2000). An artificial viscosity term was added to the pressure 
term such that 



Pi, 

„2 



1 J 



(ii) 



where the artificial viscou s pressure is given by the 
ICingold fc Monaehan I l| 19831 ) equation 



n, 



(12) 



Here we are using the notation Aij = § (Ai + Aj) and Cjj is 
the average sound speed at particle positions i, j and 



h 



(13) 



In this equation Vj 



and if 



O.Olhfi which 



prevents the denominator going to zero. We used a = 0.5 
and fJ — for the simulations reported here. 



4-1-4 The radiation force 

We follow the geometrical definitions of lPringiel (1 1996)') an d 
define a unit tilt vector (see Figure 1 of lOgilviel (|l999h ) 
1(7?, t) such that 



1 = (cos 7 sin /3, sin 7 sin /3, cos /3) 



(14) 



where /3(R,t) and ~{(R,t) are the Euler angles of the tilt 
vector with respect to a fixed Cartesian coordinate system 
(OXYZ) centered on the compact object. The position vec- 
tor of a point on the disc with distance R from the origin 
and azimuth <j> is given by Ren where eR is the radial unit 
vector and is given by 



e R = 



cos <j> sin 7 + sin <j> cos 7 cos ft 
— cos <j) cos 7 + sin <f> sin 7 cos /3 
— sin ij> cos /3 



(15) 



The radiation source is assumed to be centered on the com- 
pact object position and radiates isotropically. The radiation 
flux at a distance R from the compact object is given by 



f = 



4nR 2 

and the power absorbed at a surface element dS is 
dP = (1 - 77) If -dS\, 



(16) 



(17) 



where L* is the total luminosity of the radiation source, r/ 
is the disc albedo, that is the amount of radiation scattered 
by the surface of the disc, in all the simulations reported in 
this paper rj was assumed to be zero, and dS is an element 
of the disc surface area given by 



(18) 



We assume that the radiation is absorbed by this element 
and is uniformly re-radiated immediately and on the same 
disc side. The force due to absorption of the radiation is 



, A dP 
dA = e R , 

c 



(19) 



where c is the speed of light. 

Since the incident radiation is parallel to the radius vec- 
tor a circular disc annulus receives no torque from the ab- 
sorbed radiation. There is, however, a net torque on a circu- 
lar annulus due to the summed radiation pressure reaction, 
as calculated by iPringle I (| 19961 ). Since the radiation is re- 
radiated uniformly from the disc surface, the element will 
receive a radiation pressure reaction from this radiation of 
the form 

9 HP 

dF = -f^n, (20) 
3 c 

where n is the unit normal pointing away from the disc 
surface that received the initial radiation. 



5 SIMULATIONS 

The results presented here were generated using a local pri- 
vate Linux system and The Swinburne Centre for Astro- 
physics and Supercomputing facility. 

In the simulations the total system mass, Mt, and the 
binary separation, a, were both scaled to unity. The binary 
orbital period, P or b, was scaled to 2n. The radiation source 
was also scaled such that 



L, 



Physical q2/3 

/ 2nM t 



V Po 



5/3 



(21) 



In the simulations the accretion disc had an open in- 
ner boundary condition in the form of a hole of radius 
n = 0.0125a centred on the position of the primary ob- 
ject. Particles entering the hole were removed from the sim- 
ulation. Particles that re-entered the secondary Roche lobe 
were also removed from the simulation as were particles that 
were ejected from the disc and had a distance > 0.9a from 
the centre of mass of the system. 

We assumed a local isothermal equation of state and 
that the dissipation was radiat ed promptly from the point 
at which it was generated. The IShakura h, Sunvaevl (|l973T ) 
viscosity parameters were set to cti ow of 0.01 and cthigh of 
0. 10, and the v i scosity state changed smoothly as described 
in lTruss et"ai~l i|2000n . The SPH smoothing length, h, was 
allowed to vary in both space and time and had a maximum 
value of 0.01a. 



5.0.5 The gas stream 

We simulated the mass loss from the secondary by introduc- 
ing particles at the inner Lagrangian point (Li). See Paper 
I for a description of the particle initial velocities. 



5.0.6 The initial non-warped accretion disc 

The simulations were started with zero mass in the accre- 
tion disc and with the radiation source switched off. A sin- 
gle particle was injected into the simulation every 0.0lQ.~^ b 
at the L\ point until a quasi-steady mass equilibrium was 
reached within the disc. The simulations were then contin- 
ued for another 3 orbital periods to ensure mass equilibrium. 
The number of particles in each accretion disc was approxi- 
mately 30, 000 giving a good spatial resolution. All the sys- 
tems modelled had moderate vertical resolution of the order 
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3-5 smoothing lengths. The radiation source intensity was 
turned on gradually using a time exponential that increased 
the power of the radiation over a short period of time. 



all other input parameters were adopted from the literature 
before the simulations were run. 



5.1 Surface finding algorithm &: self-shadowing 

The accretion disc surface finding algorithm has been de- 
scribed in Paper I. Basically, the method uses a convex hull 
algorithm to find the surface particles. The accretion disc 
surface was then constructed from the set of surface points 
that form a convex hull of the disc particles. We used a ray- 
tracing algorithm to determine regions of self-shadowing. 
For each particle found on the disc surface a light-ray was 
projected from the particle to the position of the radiation 
source at the centre of the disc. The particle was deemed 
to be illuminated by the radiation source if this light-ray 
did not intersect any disc material between the particle sur- 
face position and the radiation source. The radiation force 
was only applied to particles that were considered to form 
part of the disc surface and were illuminated by the central 
radiation source. 



5.2 Disc warping and precession measure 

The two measures defined bv lLarwood fc Papaloizoul (|l997l ) 
were used to measure the disc warping and the amount of 
warp precession. They defined an angle j as the angle be- 
tween the total disc angular momentum vector and the an- 
gular momentum vector for a specific disc annulus, i.e. 

J a • Jd 



COS J 



(22) 



The term Ja is the total angular momentum within the 
specific annulus and was calculated by summing the angular 
momentum for each particle within the annulus. The term 
Jd is the total disc angular momentum and was calculated 
by summing all the angular momenta for all particles within 
the disc. An angle II was also defined which measures the 
amount of precession of the disc angular momentum relative 
to the initial binary orbital angular momentum, Jo: 



cos n = 



(Jp x Jp) ■ u 
|Jo x JdIIuI 



(23) 



where u is any arbitrary vector in the binary orbital plane. 
We also measured the tilt of the entire warped accretion disc 
relative to the orbital plane using 

Jo • Jo 



cos f3 - 



Jo Jd 



(24) 



5.3 Parameters Adopted and Adjusted 

For each simulation we adopted system parameters from the 
literature, as summarised in Tables 1, 2, 3 and 4. We built up 
the disc to mass equilibrium before turning on the radiation 
source. The precession rate in the simulations was dependent 
on the radiation intensity, as described in Paper I; by adjust- 
ing the value of L* we were able to reproduce the observed 
disc precession rates. We emphasize that the ONLY 
parameter we adjusted to match the published pre- 
cession rate was the luminosity of the central source, 



6 SIMULATION RESULTS 
6.1 HerX-1 

The SPH accretion disc model for the HerX-1 simulation 
without the radiation source was roughly point-symmetric 
about the compact object. The disc did not precess and re- 
mained approximately constant in shape and size once mass 
equilibrium had been reached. 

When the radiation source was switched on a strong 
twist developed in the disc and the entire disc tilted out of 
the orbital plane. Figure [T] contains projection plots for the 
HerX-1 model; the position of each SPH particle is indicated 
by a small black dot. The upper left hand plot, labelled xy- 
view, is a plan view of the accretion disc as seen from above 
the disc. The solid dark line is the Roche lobe of the pri- 
mary and the L\ point is to the right of the plot. Material 
from the secondary enters the primary Roche potential from 
the L\ point. The cross at the centre of the map is the posi- 
tion of the primary object. The disc remained approximately 
constant in shape and size throughout the simulation. 

The two upper right-hand plots of Figure [TJ labelled yz- 
view and xz-view, are side views of the disc in the y-z and x-z 
directions respectively. The yz-view plot is a projection view 
of the disc as seen from the secondary, similarly the xz-view 
is a projection plot with the secondary located to the right 
of the plot. The lower plot of Figure [T] labelled accretor- 
view, shows the distribution of the particles as seen from 
the compact object. The disc tilt out of the orbital plane 
is clearly seen in the yz-view plot, where the disc is viewed 
almost side on. The inner disc warp can be seen on either side 
of the position of the accretor. This warp is odd symmetrical 
about the centre of the disc. The maximum value of the warp 
is located at a distance of approximately 0.1a either side of 
the primary position. 

Figure [T] shows that the radiation force has pushed the 
entire disc out of the orbital plane. The inclination of the 
disc remained constant for many orbital periods and the disc 
precessed in a retrograde direction relative to the orbital 
flow. 

In Figure[2]we show the shape of the disc viewed at dif- 
ferent orientations. Note the z-scale, i.e. the direction per- 
pendicular to the orbital plane, is magnified by a factor of 
two, making the vertical structure more visible. The tilt an- 
gle visually inferred from Figure is therefore misleading. 
The disc is illuminated by an external light source that lies 
in the z-direction. The curved white line in each plot is the 
Roche lobe of the primary object and the straight white 
lines indicate the xyz axes defined above. The central image 
(e) shows the accretion disc as viewed from above, the disc 
shape and warp is clearly seen in this image. The images 
(b), (d), (f) and (h) are views of the disc from the sides, 
back and front of the disc. The remaining images (a), (c), 
(g) and (i) are views of the disc from slightly above or below 
the orbital plane. 

Figure [3] is similar to the xy-view in Figure [1] Here 
particles that lie above the orbital plane are indicated by a 
small square that appear dark on the plot. All other parti- 
cles are indicated using a small dot. The boundary between 
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Normalised binary separation 

-0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 -0.4 -0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Normalised azimuth angle (radians/2it) 



Figure 1. Particle projection plots for HerX-1 simulation. The position of each particle is indicated by a small black dot. The plot 
labelled xy-view is a plan view of the accretion disc as seen from above the disc. The cross at the centre of the plot shows the position 
of the primary object. The solid dark line is the Roche lobe of the primary and the L\ point is to the right of the plot. The two plots 
xz-view and yz-view are particle projection plots on a plane perpendicular to the orbital plane and through the system axis. The bottom 
plot, accretor-view, shows the particle distribution as seen from HerX-1. The horizontal axis is the normalised azimuth angle, the L\ 
point is at azimuth angle and the stream/disc impact region is at an angle of approximately 0.9. The vertical axis is the angle, in 
degrees, between a particle and the orbital plane when viewed from the compact object. The disc material flows from right to left. 



the two sets of points shows the Eulerian angle of the warp 
precession, 7. For this snapshot the Eulerian angle, 7, is 
approximately odd symmetrical with respect to the y-axis. 
This boundary line describes the intersection between the 
warped accretion disc and the orbital plane. As the disc pre- 
cesses in a retrograde direction the position of the Eulerian 
angle, 7, also rotates at the same rate. 

The disc precessed in a retrograde direction relative to 
the inertial frame, the entire disc precessing with a single pe- 
riod. This precession is accompanied by continuous flexing of 
the disc structure as its orientation with respect to the Roche 
potential changes due to the orbital motion of the donor star. 
We tuned the luminosity to a value of L* = 2.15x 10 37 erg s _1 
to achieve the observed precession rate of approximately 34.9 
days. This value of L» is approximately 10% higher than the 
value quoted in Table 2. This discrepancy is, of course, well 
within the uncertainty arising from the distance adopted by 



iMcCrav et al~1 l|l982l ). More importantly, the factor implic- 
itly introduced by our assumption of zero albedo is proba- 
bly approximately 2, i.e. far larger than approximately 10%. 
The warp shape was very stable over many orbital periods; 
we ran this simulation for more than 60 orbital periods and 
the warp shape remained constant throughout. The rate of 
precession of the disc was determined as described in paper 
I. Figure U shows the precession angle versus orbital phase 
for the model. The disc precession rate was determined by 
fitting a straight line t o the data using a Numerical Recipes 
least squares method (|Press et al. Ill986l ). 

The warped disc is inclined relative to the orbital plane 
as shown in Figure [T] (yz-view) which changes the point at 
which the gas stream impacts the accretion disc. For a planar 
accretion disc the material from the gas stream always im- 
pacts the outer edge of the disc at approximately the same 
point and the energy dissipated at this point is approxi- 
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Figure 2. Her X-1 SPH accretion disc viewed at different angles. Note the direction perpendicular to the orbital plane has been artificially 
magnified by a factor of two. This renders the vertical structure more visible, but causes the tilt angles to appear much larger than they 
really are. The curved white line is the Roche lobe of the compact object and the straight white lines form axes centred on the compact 
object. The disc is illuminated using an external white light source place along the z-axis. 



mately constant with time. This is not true for a warped 
disc. As the secondary orbits around the warped, tilted disc, 
the impact point of the gas stream moves from the edge of 
the disc to nearer its centre. Consequently the kinetic energy 
gained by the ballistic stream, and subsequently dissipated 
at the impact, varies. 

The lower plot of Figure [5] shows how the distance 

from the compact object of the stream/disc impact point 

varies with the orbital phase. This plot does not repeat with 

the orbital phase, instead it repeats with a period of 1.65 

days which is the synodic period between the orbital period 

and the retrograde precession period of the accretion disc: 
p-i _ p-i I p-i 

± syn ± or l) I ± prec- 

The upper plots of Figure [5] (a, b, c and d) show disc 
density maps and the ballistic trajectory of the gas stream 
from the L\ point to the impact point on the accretion disc. 
The orbital phases of plots (a, b, c and d) are indicated on 
the lower plot. Plots (a) and (c) show that the gas stream 
impacts the disc very close to the central object. The density 
maps also show that the disc is eccentric and changing shape 
on the orbital period. The tidal action of the secondary 



is making the disc flex and distort. Tidally-driven flexing 
will occur whenever the disc is no t fixed wit hin the pri- 
mary Roche lobe; see lHaswell et al~l l|200lf ) and lSmith et al.l 
(2007) for a discussion of this for apsidally precessing discs. 



In the interval between orbital phases approximately 
0.11 and approximately 0.31 the gas stream flows over the 
disc surface and impacts the disc very close to the central 
object (see Figure [5]). At orbital phase approximately 0.31 
the impact point very abruptly changes to near the disc edge. 
The stream continues to interact with the edge of the disc 
for approximately 0.18P o rt- The gas stream then interacts 
with the other surface of the disc and slowly the warp and 
tilt opens allowing the gas stream to flow further towards 
the central object. There is then a second abrupt change of 
the impact distance at about orbital phase approximately 
0.75 as the gas stream again impacts close to the edge of 
the disc. This change is not as large as the first because 
the disc shape is asymmetric. The gas stream interacts with 
the disc near the edge of the disc for approximately 0.16 
Porb and then slowly moves closer to the central object on 
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Figure 5. Penetration of the gas stream for the HerX-1 simulation. The lower plot shows the stream/disc impact distance from the 
neutron star as function of orbital phase. The four upper plots are density maps of the accretion disc as viewed from above of the disc. 
The Roche lobe of the neutron star is also plotted. The ballistic path of the gas stream is plotted, it starts at the L\ point and ends 
when the gas stream intercepts the accretion disc. 



the first side of the disc. The sequence then repeats on the 
synodic period. 

The simulated disc dissipation light curves for the 
warped disc are shown in Figure [6] Dissipation curve (a) 
is the dissipation rate integrated over the entire disc and 
was evaluated by summing the dissipation rate from all the 
particles in the disc. Curves (b) , (c) and (d) show the dissipa- 



tion from the disc at distances greater than 0.09a, 0.16a and 
0.25a respectively. The peak at orbital phase approximately 
0.3 corresponds to when the gas stream intercepts the accre- 
tion disc close to the compact object, labelled (a) in Figure 
this recurs on the 1.65 day synodic period, which is a few 
percent shorter than the orbital period. Most of the variabil- 
ity in Figure [S] is caused by the stream impact moving close 
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Intersection of Her X-1 disc and orbital plane 



Disc dissipation light curve Her X-1 
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Figure 3. Her X-1 particle plot plan view. The particles above 
the orbital plane are indicated using a small square that appear 
dark in the plot. The remaining particles, that are below or on 
the orbital plane are indicated using a small dot. The boundary 
of the two sets of points shows the intersection of the accretion 
disc and the orbital plane. 



Precession angle versus orbit phase for Her X-1 




Orbital phase 



Figure 4. Disc precession angle versus orbital phase for the 
Her X-1 SPH simulation. The graph shows the direction of the 
angular momentum for the whole disc measured relative to some 
arbitrary starting position. The disc rotates in a retrograde di- 
rection with a precession period of approximately 34.9 days. 



to the compact object. Figure 6 has no direct observational 
counterpart: for the Her X-1 system we expect the observed 
optical light curve to be dominated by reprocessing of X- 
rays generated very close to the compact object; modelling 
this reprocessing-powered optical light curve is beyond the 
scope of this paper as it would require realistic treatment of 
the wavelength dependent albedo. 

A radiation source centred on the compact object will 
be hidden to a high inclination observer as the tilted disc 
moves between the observer and the radiation source. When 




Orbital phase 

Figure 6. Her X-1 dissipation rate light curves for different re- 
gions of the accretion disc. Curve (a) is the dissipation rate in- 
tegrated over the entire disc. Curves (b), (c) and (d) show the 
dissipation from the disc at distances greater 0.09a, 0.16a and 
0.25a respectively. 



the line of sight is clear the central X-ray source is visible, re- 
sulting in an X-ray on state. However, when the disc material 
comes between the observer and the X-ray source the X-rays 
will be absorbed by the disc material leading to an X-ray off 
state. In Figure [7] we show curves marking X-ray on and off 
states as seen from different observer inclinations. For the 
first three panels where the inclination is less than 82° only 
one on and off state occurs per precession period. For larger 
inclinations there are two on states of different lengths, one 
when the line of sight passes over the top of the disc and an- 
other when the line of sight passes under the disc resulting in 
a short on state. X-ray observations of Her X-1 show similar 



behaviour (IDeeter et al. 



198ll; Middlcdi tch fc Nels"on~lll976l ; 



iGerend fc Bovnton II 19TE ) with a "main-high" and a "short- 
high" on states. The inclination of i > 82° implied by our 
model is consistent with independen t estimates from obser- 
vations of Her X-1. Recent pape rs e.g. Boroson. et al. I (I2007T) 
state simply i > 80° , and cite iHowarth fc Wilson! (|l983l Vs 
value of i = 83.25° ± 0.25° though the latter state "we have 
reservations about the quoted inclination and disc tilt" and 
"it may have been preferable to adopt a fixed value for the 
inclination (i ~ 85°)". Though we recognise the limitations 
of our model, we are encouraged that by tuning a single 
parameter, Lx, to reproduce the disc precession period, we 
simultaneously reproduce the observed X-ray luminosity, the 
X-ray on/off behaviour and a consensus value of the orbital 
inclination. We conclude that our work is strong evidence 
that irradiation-driven warping is the physical cause under- 
lying the precessing tilted disc in Her X-1. 



6.2 SS433 

Using parameters from Table [2] the accretion disc for the 
SS433 model without the radiation source was also roughly 
point-symmetric about the compact object. The disc did not 
precess and remained approximately constant in shape and 
size once mass equilibrium had been reached. When the ra- 
diation source was switched on a strong twist developed in 
the disc and the entire disc tilted out of the orbital plane. 
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Figure 8. Particle positions plots for SS433. See caption of Figure[T] 



Figure [8j shows particle projections (as described in Figure 
[1} for the SS 433 simulation. The plan view of the accretion 
disc shows the disc to be asymmetric. The disc had two spi- 
ral density compression waves at the edge. The plan view of 
the disc remained approximately constant throughout the 
run. The two upper right-hand plots of Figure [8j show that 
this disc was inclined with respect to the orbital plane. In 
the yz-view the disc is viewed almost side on. This leads 
to the gas stream impacting the disc near the circulariza- 
tion radius instead of the edge of the disc. The inner disc 
regions are relatively close to the orbital plane, while the 
out er parts of the disc are ti lting to large angles (c.f. Figure 
7 of lWiiers fe Pringle |[l999l ). The lower plot of Figure[8l la- 
belled accretor-view, shows the distribution of the particles 
as seen from the accretor. As in the case of Her X-l, the radi- 
ation force has tilted the entire disc out of the orbital plane. 
The inclination of the disc remained constant for many or- 
bital periods and the disc precessed in a retrograde direction, 
see Figure [9] We ran this simulation for 50 orbital periods 
in which the tilt and warp shape remained approximately 
constant. With a central luminosity of 1.2 x 10 36 ergs -1 , we 
obtained a precession rate of approximately 162 days. This 
luminosity is a factor of approximately 4 above the 2-10 keV 



luminosity found bv lMarshall et al. I (I2002T) and about twice 
the 2-8 keV luminosity found bv lKotani et al.~l (|l996h : both 
these observed values are sensitive to the spectral model 
and -/V H . In Table [2] we quote the (model-dependen t ) lum i- 
nosity from iBrinkmann. Kawai. Matsuoka fc Fink I (|l99lf ). 
which is 20% less than the Lx we obtain by tuning to the 
disc precession period. More important than these discrep- 
ancies in luminosity is the kinetic power in the jets, which 
is approximately 1000 times the X-ray luminosity. We sur- 
mise that this will have a significant effect on the dynamics 
of the inner disc, but any detailed analysis requires a more 
definite picture of the jet acceleration mechanism than has 
yet emerged. As noted earlier in our discussion of Her X-l, 
our assumptions about the albedo must also be considered 
when comparing our model luminosities to the observed lu- 
minosities: since we assume zero albedo, we are effectively 
boosting the irradiation by a factor of approximately 2 at a 
given luminosity. 

This warped disc is also inclined relative to the orbital 
plane as shown in Figures [51 and fTOl As illustrated for HerX- 
1 in Figure O as the secondary in SS433 orbits the tilted 
disc, the impact point of the gas stream varies from the edge 
of the disc to nearer its centre. 
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Figure 10. SS433 SPH accretion disc viewed at different angles, the z values have been artificially magnified by a factor of two, so the 
tilt angles here are misleading. See caption of Figure [2] 



The simulated disc dissipation light curves for the 
SS 433 warped disc are shown in Figure[TT] Dissipation curve 
(a) is the dissipation rate integrated over the entire disc. 
Curves (b), (c) and (d) show the dissipation from the disc 
at distances greater 0.09a, 0.16a and 0.25a respectively. The 
dip at orbital phase approximately 0.55 corresponds to when 
the gas stream intercepts the accretion disc edge. The peaks 
at orbital phase approximately 0.4 and approximately 0.7 
are due to gas stream impacts close to the central object. 
The peaks and dips repeat on the synodic period, as they 
did for HerX-1, and again most of the variability in dissipa- 
tion is generated by the stream impact moving close to the 
compact object. 

The overall disc tilt out of the orbital plane was found 
to be 11.3° with a maximum warp tilt angle of 18.6° near the 
edge of the disc. The tilt of the disc will lead to the radiation 
source being hidden to a high inclination observer as the 
tilted disc moves between the observer and the radiation 
source. This will occur as the warped disc precesses, as it 
did in the HerX-1 simulation (see Figure [7J. The radiation 
source may also be obscured by the material in the warp 
close to the centre of the disc. 

By tuning the central luminosity we have been able 



to reproduce the observed long period for SS433. The tilt 
angle of the disc at the edge is very similar to the half- 
angle (20°) for the collimated precessing jets. In the cen- 
tre of the disc, which is presumably where the jets are 
launched, the disc tilt angle was, however, about 12°. It is 
intriguing to wonder whether this discrepancy might be di- 
minished with the system parameter s recen tly published by 
iBlundell. Bowler, fe Schmidtobreickl l|2008h . We note also, 
however, that there is evidence for an accretion disc wind 
from this system, which along with the acceleration of the 
relativistic jets, is absent in our model. 



6.3 Other systems modelled 

We also modelled eight other systems, Tables [3] and U con- 
tain the main system parameters adopted in this study for 
these systems. The eight models simulate five X-ray bina- 
ries with super-orbital periods, th e source 4U 1626-67 dis- 
cussed bv lWiiers &: Pring lc (1999), a 'generic' LMXB, and 
the CV KRAur. For the 'generic' LMXB, Pi ong is that of 
our model. KRAur i s a 'negative superhump' system which 
iKozhevnikoTI |2007l ) suggests may harbour a radiation- 
driven warped disc. In each relevant case, the model lumi- 
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Figure 7. X-ray on and off states for two complete precession 
periods viewed at different inclinations for the HerX-1 simulation. 
In the graphs a value of 1 indicates an X-ray on state, a value of 
is an X-ray off state. The observation angle is indicated in the 
upper left of each plot. 



Precession angle versus orbit phase for SS 433 




Orbital phase 

Figure 9. Disc precession angle versus orbital phase for the 
SS 433 SPH simulation. The graph shows the direction of the an- 
gular momentum for the whole disc measured relative to some 
arbitrary starting position. The disc rotates in a retrograde di- 
rection with a precession period of approximately 162 days. 



nosities were chosen to reproduce the observed long periods 
found in these systems. In all cases the disc developed strong 
warps, detailed in Table [5] 

For systems with extreme mass ratios, i.e. where 
M C ompact > approximately 4Mdonor (c.f. ISmith et al.ll2007h 
a strong warp develops close to the centre of the disc to the 
extent that some of the outer parts are completely shadowed 
by the inner warp. The warp precesses in a wave-like man- 
ner in a retrograde direction relative to the inertial frame. 
In addition, the outer disc precesses in a prograde direc- 
tion due to the well-known eccentri c instability first discov- 
ered in the SUUMa CVs (see e.g. IVogtJll982i : ISmith et all 




Orbital phase 



Figure 11. SS433 dissipation rate light curves for different re- 
gions of the accretion disc. See caption of Figure [6] 



l2007ll In our extreme mass-ratio warped disc simulations 
(i.e. those for X 1916-053 and 4U 1626-67) superhumps are 
apparent in the simulated dissipation light curve, just as 
they were for the plane disc simulated by iFoulkes et al. I 
(2004). We found a similar behaviour in the simulations of 
Bo 158 l|Barnard et al hood ). For the less extreme mass ratio 
systems a small warp is induced in the inner regions of the 
disc. This warp does not shadow the outer regions of the disc 
and these regions also warp. The entire disc then becomes 
tilted to the binary plane and the whole disc precesses in 
a retrograde direction. This is reflected in Table [5] the two 
extreme mass ratio systems have their maximum warps in 
the inner disc. Generally the other systems have their max- 
imum warp at the disc edge, but the CygX-1 model is an 
exception to this rule. 

4U 1626-67 was modelled to examine the hy pothesis of 
iKerkwiik et al. I i|l998t ). IWiiers fc Pringiel l| 19991 ) and others 
that for large central luminosities the inner disc warp angle 
exceeds 90° , so these disc annuli flow in a direction counter 
to the orbital angular velocity. If this were the case, the 
reversal of the disc angular momentum vector could account 
for th e torque reversals observed by e.g. IChakrabartv et al. I 
(|l997h . As Table \5\ shows, we find an inner disc warp of 21° 
for 4U 1626-67. This is the largest angular deviation from the 
orbi tal plane found in any o f our simulations, so we agree 
with lWiiers fc Pringle I (|l999h that 4U 1626-67 is an extreme 
case. We suggest, however, that the inversion of the disc 
angular momentum vector in their inner disc is simply an 
artefact of the approximations necessitated by their analytic 
approach. Our conclusion leaves the mechanism underlying 
the observed torque reversals as an open question. 

The luminosity used in each system model was fine- 
tuned to reproduce the observed long periods and tabulated 
in Table [5] For all eight X-ray binaries with observed super- 
orbital periods, the luminosities required to produced the 
observed long periods are similar to the observed luminos- 
itjEI, see Figure [12] The biggest discrepancies are about 20% 

1 In the tables and Fig. 1121 we have used the long-term average 
luminosity, derived from the RXTE ASM, where this has been 
published. For three of the four brightest sources the error on 
this value is smaller than the size of the symbol used. For the 
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Table 3. Major parameters for the other binary systems modelled 



System 


^orb 


q 


Aii 


M s 


-flong 


Reference 


name 


(days) 


M 2 /Mi 


M 


(lO-SMeyr- 1 ) 


(days) 




SMC X-l 


3.9 


12.9 


1.4 


4.0 


« 55 


1 


Cyg X-l 


5.6 


2.1 


16.0 


2.5 


294 


1 


Cyg X-2 


9.844 


0.5 


1.4 


1.0 


78 


1 


X 1916-053 


0.035 


0.07 


1.4 


0.1 


3.8 


1 


4U 1626-67 


0.029 


0.02 


1.4 


0.05 


0.0 


1 


LMC X-3 


1.7 


0.6 


10.0 


3.0 


198 


1 


KR Aur 


0.16 


0.6857 


0.7 


1.0 x 10 -10 


5.3 


2 


Gen LMXB 


0.2 


0.36 


1.4 


0.1 


5.25 


1 



Reference: 1 IWiiers fc Pringle I d 199911 , 2 iKozhevnikoTl J2007h 



Table 4. X-ray Luminosity and system type 



System 


X-ray Luminosity L* 


LEdd 


Type 


Reference 


name 


(xlO 38 ergs" 1 ) 


(x 10 38 ergs -1 ) 






SMC X-l 


2.0 


1.82 


HMXB 


1 


Cyg X-l 


0.05 


20.8 


HMXB 


1 


Cyg X-2 


1.8 


1.82 


LMXB 


1 


X 1916-053 


0.05 


1.82 


LMXB 


2 


4U 1626-67 


0.1 


1.82 


LMXB 


3 


LMC X-3 


1.5 


13.0 


HMXB 


1 


KR Aur 


see text 


see text 


LMXB 


1 


Gen LMXB 


1.0 


1.82 


LMXB 


1 



Reference: 1 iGrimm. Gilfanov fc Sunvaev I J2002I ), 2 lNarita et aTTI (l2003h ■ 3 IChakrabartv:1998 I lll998h . 4 iKozhevnikov I ll2007l) 



in the cases of SS433 and SMC X-l. As already noted, our 
model is of limited applicability to SS433 because we do 
not attempt to model the kinematic luminosity. SMC X-l is 
probably partly powered by wind-fed accretion, altering the 
angular momentum balance in the disc from that in our ex- 
clusively Roche lobe overflow fed model disc. Given the sim- 
plistic treatment of the reprocessing and the limitations in 
our simulation of the accretion disc viscosity, this agreement 
is astonishingly good. We conclude that irradiation driven 
warping is probably generally the mechanism underlying the 
observed long periods in X-ray binaries. 

For KR A ur the X-ray luminosity is poorly known. 
IShafterl (|l983l ) cites a private communication giving the Ein- 



This i mplies a L x 



erg sec 1 using a distance of 



180 pc (|Greiner Ill998l ). The ASCA 0.7 - lOkeV blackbody 
fit (a poor fit to the spectrum) gives Lx = 10 30 ' 4 erg sec -1 
using a distance of 100 pc or Lx = 10 30 ' 9 erg sec -1 using 
a distance of 180 pc. These numbers are very low com- 
pa red to the lum i nositi es mentioned in the final discussion 
m IKozhevnikoTl (|2007l ). and are indeed typical of the X- 
ray /optical flux ratio s of optically selected CVs (|Shafterl 
1 19831 : Patterson Ill98ll ) . 

As Table [5] shows, the model Lx required to match 
the inferred warp disc precession period in KRAur is 



fourth, we used the data of Staubert et al (2007) to deduce the 
size of the error bar, again finding it to be smaller than the symbol 
used. For the three least luminous sources we were unable to a 
reliable error bar for the long-term average luminosity. 



10 37 erg sec 1 . Intriguingly, thi s is typical of the X-ray lu- 
minosity for supersoft sources remcr l l2000l ). 



6.4 Orbital inclinations and X-ray Dips 

In LMXBs X-ray dips are well-understood phenomena, aris- 
ing when material from the accretion flow transiently blocks 
our line of sight to the central X-ray source. Table[S]gives the 
maximum warp angle and overall disc tilt angle for each sys- 
tem simulated. These results allows us to place constraints 
on the systems' orbital inclinations using observations of 
the presence or absence of X-ray dips. Or, alternatively we 
can use the presence or absence of X-ray dips, coupled with 
known orbital inclinations, to examine whether our results 
make sense. Her X-l and SMC X-l are eclipsing, so the in- 
clinations are strongly constrained by this. We have already 
shown that the disc shape we find for Her X-l is consistent 
with the X-ray on and off behaviour for i « 82° . 

Of the LMXBs we treat, all show X-ray dips. Since we 
selected them for their observed super-orbital periods and 
well-constrained system parameters (except for 4U 1626-67, 
selected because previous authors discussed it in the con- 
text of warped precessing discs) this is interesting. It implies 
that the long period modulations may be preferentially de- 
tected at high orbital inclinations. This is consistent with 
the moderate values we find for the angular deviation of the 
disc from the orbital plane. The biggest overall disc tilt is 
12.3° for SMC X-l while the most pronounced angular de- 
viation is the 21° inner warp in 4U 1626-67 (see Table [5]). 
At low inclinations the disc will always appear more-or-less 
the same if the overall disc tilt has the moderate values we 
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Table 5. Results the other binary systems modelled. Columns are: 1 System modelled, 2 Model luminosity used to reproduce the system 
long period, 3 Maximum disc warp angle, 4 Position of the maximum disc warp, Inner disc or Edge of the disc and 5 The overall disc 
tilt. 



System 
name 


Model Luminosity L* 
(x 10 38 ergs^ 1 ) 


Maximum Disc Warp 
(Degrees) 


Maximum Disc 
Warp Position 


Disc tilt 
(Degrees) 


SMC X-l 


1.65 


20.9 


Edge 


12.3 


Cyg X-l 


0.05 


14.7 


Inner 


6.2 


Cyg X-2 


1.70 


19.4 


Edge 


10.5 


X 1916-053 


0.05 


12.8 


Inner 


0.8 


4U 1626-67 


0.1 


21.0 


Inner 


1.0 


LMC X-3 


1.75 


19.6 


Edge 


10.6 


KR Aur 


0.1 


7.5 


Edge 


3.0 


Gen LMXB 


1.0 


18.5 


Edge 


10.6 


Her X-l 


0.215 


13.7 


Edge 


5.5 


SS433 


0.012 


18.6 


Edge 


11.3 



find. In contrast, some of the models discussed in the lit- 
erature for Her X-l have tilts of 40°: if this were typical, 
then the precession would produce observable effects at low 
and moderate inclinations. As a caveat we note that system 
parameters are generally easier to determine in high inclina- 
tion systems, so the preponderance of dippers may be partly 
due to this selection criterion. 

4U 1916-053 is an extensively-studied dip source, with- 
out a well-known orbital inclination. Our results, with a 
maximum disc warp of 12.8°, suggest an orbital inclination 
of approximately 77°, such that the inner disc warp inter- 
mittently blocks the line of sight to the central X-ray source. 
4U 1626-67 is similarly a well-known dip source. In this case, 
we find a pronounced inner disc warp of 21°, the highest for 
all our simulations. This implies that the orbital inclination 
in this case is less than approximately 70°. 

Cyg X-2, conversely, is not a well-known X -ray dip 
source. Dips were reported bv lVrtilek et al. I (|l986l ), but de- 
spite Cyg X-2 being one of the brightest persistent LMXBs 
we did not find other dips in the literature. This suggests 
Cyg X-2 has an inclination which is just lower than that 
required for dipping to be seen, and only when stochastic 
variations happen to conspire to move material an unusual 
distance above the disc mid-pla ne do dips occur. The incl i- 
nation of Cyg X-2 was studied bv lOrosz fc Kuulkers ll| 19991) . 
They find i ~ 55°, from modelling of the ellipsoidal varia- 
tions assuming the disc has a temperature profile T <x r _3//r 
typical of steady irradiated discs; or i = 62.5° ±4° assuming 
the canonical steady-state disc T oc r _3//4 . Coupled with our 
maximum disc tilt of 19.4°, these numbers are consistent 
with dipping only in rare, extremely favourable stochastic 
deviations. 

Taken overall, therefore, we consider our results for the 
LMXBs to be encouragingly consistent with the observa- 
tional literature. Emboldened by this, we take our values 
for the maximum disc warps in 4U 1916-053 and 4U 1626-67 
at face value, and venture to predict that the inclinations 
in these two systems are approximately 77° and approxi- 
mately 69°. These are simple estimates, ignoring the effect 
of the finite opening angle of the disc. If the warped disk 
maintains the canonical H/R ~ 0.1 relationship, then our 
inclination estimates would be reduced by approximately 6° 
in each case. 



We also examined the literature for dips in the HMXBs, 
finding nothing except in the case of Cyg X-l. In this case the 
orbital inclination is a ccepted to be approximately 30° (e.g. 
|Gies fc Bolton 1 1 1986] )) and the dips are likely associated 
with material in the wind from the donor star. Where discs 
are wind-fed, the specific angular momentum and the X-ray 
luminosity is generally less. Thus HMXB discs are generally 
smaller and less irradiated than LMXB discs, both charac- 
teristics tending to reduce the angular elevation of outlying 
material above the disc mid-plane. Our simulations show 
that for HMXB parameters, the discs invariably have their 
maximum warps at the outer edge, lacking the pronounced 
inner disc warps of the extreme mass ratio LMXBs. All of 
these factors conspire to make dips from accretion disc ma- 
terial in HMXBs less prevalent than in LMXBs. 



7 CONCLUSIONS 

We have studied the radiation-driven warping of thin accre- 
tion discs in X-ray binaries and one CV using a non-linear 
SPH code. A convex hull algorithm was used to find the sur- 
face particles of the disc and a ray-tracing algorithm used to 
evaluate regions of self-shadowing, as described fully in Pa- 
per I. Initially planar accretion discs were illuminated with 
radiation from the centre of the discs. The surface of each 
disc was subject to a force from the absorption of the X-rays 
on the disc surface and back reaction from the remitted ra- 
diation. 

The radiation caused these discs to warp and tilt out of 
the orbital plane, precessing in a retrograde direction rela- 
tive to the inertial frame. The resulting tilt and precession 
rate remained constant for many orbital periods. 

We simulated a number of systems that have long peri- 
ods in their observed light curves. We tuned the irradiating 
luminosity, L», to match the precession-rate of the warped 
disc to the observed long periods in seven X-ray binaries 
and the inferred warped disc precession period in one CV. 
Figure [12] shows the astonishingly good match between ob- 
served Lx and the irradiating luminosity, L, , demanded by 
our models to match the observed periods. The largest de- 
viations from the 1:1 straight line in Figure [12] are approxi- 
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mately 20%, less than the typical uncertainty introduced by 
the spectral model and the d istance determination. 

lOgilvie fc Dubus I (l200lh [OD01] found that only low- 
mass X-ray binaries with high orbital periods are likely to be 
unstable to irradiation-driven warping. In contrast to this, 
our simulations find irradiation-driven warping for all the 
orbital periods considered. In common with us, OD01 con- 
sidered self-shadowing and non-linear solutions. Their work 
is analytical and consequently necessarily makes approxi- 
mations. Our simulations are able to capture more complex 
behaviour. For example OD01 state "Owing to the symme- 
tries of the problem, these non-linear solutions are steady in 
the rotating frame of reference..." whereas our disc contin- 
uously flexes in response to the changing orientation of the 
Roche potential. Our short orbital period LMXBs tend to 
have extreme mass ratios, and for these the motion is par- 
ticularly complex, with prograde apsidal precession in the 
outer disc, warped retrograde precession in the inner disc 
and the spiral structure all contributing to a changing self- 
shadowing of the disc surface. Using analytical techniques it 
is unfeasible to robustly capture these complex motions on 
a variety of timescales. 

In Figures [13] and [14] we show how the disc tilt and 
the disc precession periods depend on the model luminosity, 
L* . The scatter in these two plots demonstrate that the disc 
behaviour is influenced by parameters other than luminosity 
in a complex way. This strengthens our confidence in our 
conclusion that Figure [T2] provides strong evidence in favour 
of the hypothesis that the observed long periods are caused 
by the precession of irradiation-warped discs. 

HerX-1 is probably the best-characterised of the ob- 
served systems. Our simulation of its disc appears to si- 
multaneously match several of the system characteristics by 
tuning a single parameter, the X-ray luminosity. Encourag- 
ingly, the tuned value for this parameter is consistent with 
the observed value. Our HerX-1 disc produces the observed 
"main-high" , "short-high" and off states for an orbital incli- 
nation consist ent with independent d eterminations. As re- 
cently as 2007, iBoroson. et al. I (|2007T ) stated that the phys- 
ical cause for the warped precessing disc remains unknown. 
We strongly conclude that irradiation-driven warping un- 
derlies the long-known and well-observed tilted precessing 
disc in HerX-1. 

For the other X-ray binaries, many of the parameters 
are subject to significant uncertainty. Encouraged by our 
HerX-1 results, we surmise that we may be able to use the 
results in Table 5 to constrain the system orbital inclina- 
tions, using the observed presence, absence and periodicities 
in X-ray dipping behaviour. We discussed this in section lrT4l 
finding our results to be consistent with the observations 
for all the LMXBs. We venture to predict inclinations for 
X 1916-053 and 4U 1626-67 of approximately 77° and 69° 
respectively. We also noted that we find maximum warps 
located in the inner discs for systems like these two, which 
have mass ratio extreme enough for the disc to harbour the 
3:1 resonance, leading to an apsidally precessing eccentric 
outer disc. The general lack of these pronounced inner disc 
warps for less extreme mass ratios may partially explain the 
lack of X-ray dips in HMXB systems, as discussed at the 
end of section 16.41 

We included a sim ulation of the C V KR Aur, motivated 
by the statements in iKozhevnikovl (|2007f l. We found we 




Figure 12. The observed X-ray luminosity versus the model 
X-ray luminosity required to match the observed super-orbital 
period for the seven systems with observed super-orbital periods 
(crosses). The line is Model Lx = Observed Lx- In all cases the 
model X-ray luminosity required to match the observed super- 
orbital period is extremely close to the observed X-ray luminosity. 
We consider this excellent agreement to be strong evidence in 
favour of the hypothesis that these super-orbital periods are the 
result of precession of irradiation-driven warped discs. 
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Figure 13. The model overall disc tilt versus the model X-ray 
luminosity. The scatter of points indicates that parameters other 
than Lx play an important role in determining the disc tilt. 



could reproduce the inferred tilted disc precession period 
for this system, but only by using Lx = 10 37 erg sec -1 . 
This value is about a million times higher than the observed 
values of Lx we were able to find for this system, but is 
consistent with the X-ray luminosities of supersoft sources. 
Thus our findings seem consist ent with the hypothesis ad- 
vanced by iKozhevnikov I (|2007h that disc precession could 
be present in KRAur and other VY Scl stars if they har- 
bour white dwarfs with steady nuclear burning on their sur- 
faces. The luminosity thus generated is sufficient to produce 
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Figure 14. The model disc precession period versus the model 
X-ray luminosity. The scatter of points indicates that parameters 
other than Lx play an important role in determining the disc 
precession period. 



irradiation-driven warping and tilting of the accretion disc 
which precesses and leads to the observed negative super- 
humps. The obvious test of this hypothesis is to verify that 
the inferred nuclear-burning luminosity does indeed exist. 

For HerX-1 we showed that the position of the stream 
impact on the warped and tilted disc was a function of both 
orbital and disc precession phase. Generally the gas stream 
misses the edge of the disc, which is tilted out of the orbital 
plane, as illustrated in Figures 1 and 2. The stream flows 
deeper into the gravitational well of the primary, and hits 
the face of the disc at a continually-varying distance from the 
compact object, as illustrated in Figure 5. We showed the 
resulting simulated disc dissipation light curve in Figure 6, 
but we note that for HerX-1 itself, this curve is not directly 
comparable to any observable. This is because HerX-1 is a 
relatively large binary, and the observed optical light will 
be dominated by reprocessing of the X-rays generated well 
within the closest approach of the stream-disc impact point. 

There is however, an interesting possible application of 
the varying stream impact depth: this impact may power 
X-rays directly in some very compact X-ray binary systems. 
This would require disc warping to allow the stream to pen- 
etrate close to the compact object, and is therefore similar 
in flavour to the direct-impact acc retor model discussed by 
iMarsh. Neleman fe Steeghs] l|2004 ). 
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